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We analyze hadron as well as quark and diquark correlation functions in Landau gauge in order to extract 
information on the spin dependence of the quark-quark interaction. We find evidence that the N — A mass 
splitting can be attributed to the spin dependence of the interaction between quarks in a colour anti-triplet state 
with spin and 1, respectively. The lightest excitations are observed in the S — channel. However, no evidence 
for a deeply bound diquark state is found. 
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1. INTRODUCTION 

The spin dependent interaction among quarks 
found in perturbation theory from one gluon ex- 
change diagrams (OGE) (lj as well as induced by 
instantons leads at least to a qualitative under- 
standing of the fine structure of the hadron spec- 
trum 0. In particular, the attractive interaction 
in the S = channel has led to the idea that 
diquarks may be well defined, localized objects 
(bound states?) within a nucleus. This also led to 
the speculation that diquarks may play an impor- 
tant role for the phase structure of QCD at high 
density. The possibility of a diquark gas || has 
been discussed and recently new ideas about the 
possible existence of a diquark condensate (colour 
superconductor) have been put forward. 

Whether these newly proposed phases of dense 
matter are realized in nature depends on details 
of the spin dependent part of the q-q interaction. 
Also in the case of the hadron spectrum it is of 
interest to analyze the relative importance of the 
different mechanisms that can give rise to the ob- 
served spin splitting of hadronic states. Already 
in the case of the N — A mass splitting the con- 
tributions from S = and 1 channels differ in 
the case of OGE and instanton induced interac- 
tions. Both lead to attractive terms for spin 
diquarks in a colour anti-triplet (3 C ) state they 
predict, however, opposite signs for the interac- 
tion in the S = 1 channel (see Table 1). 

While calculations within the instanton liquid 
model gave indications for a quite light (deeply 
bound) diquark state |7]] an analysis of the quark 
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Table 1 

Diquark states with spin S and different flavour 
(F) and colour (C) representations. The last 
two columns give the relative strength of inter- 
action terms corresponding to a flavour-spin (in- 
stanton) and colour-spin (OGE) coupling, Vs ~ 
(A" A2XS1S2), where A" denote SU(3) generators. 



distribution inside a nucleus did not give any 
hints for well localized diquarks || . In the follow- 
ing we will present some results from our analysis 
of the q-q interaction performed in Landau gauge 
within the quenched approximation of QCD B. 

The calculations we are reporting here are 
based on an analysis of 73 gauge field configu- 
rations, generated on 16 3 x 32 lattices at (3 = 4.1 
with a tree-level Symanzik improved action. Af- 
ter fixing the Landau gauge fermionic correlation 
functions have been calculated with four different 
source vectors using a tree-level improved clover 
action. Propagators have been analyzed for 8 val- 
ues of the hopping parameter, k € [0.14,0.148], 
which corresponds to the hadron mass interval 
0.5 < m^jmp < 0.9. From a calculation of the 
string tension we find the lattice cut-off, a -1 ~ 
1.1 GeV; our largest k- value corresponds to a pion 
mass = (0.316 ± 0.003)a" 1 ~ 350 MeV. 
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Figure 1. The N — A mass splitting versus the 
pion-rho mass ratio. The star denotes the ex- 
perimental value. Data for Wilson fermions are 
taken from [10]. Also shown is the mass splitting 
between S = and 1 diquarks (open squares). 
Errors are similar to those shown for the N — A 
mass splitting. 

2. N - A MASS SPLITTING 

On our data set we have calculated correlation 
functions for the nucleon and delta. We find that 
the N — A mass difference increases with decreas- 
ing quark massi. As shown in Fig. 1 results are 
also quantitatively consistent with earlier findings 
from calculations with Wilson fermions [lC| ]. 

Simple constituent quark models relate the N— 
A mass splitting to the spin dependence of the 
interaction among quarks in a 3 C state, i. e. iba — 
tojv = 2(Vs=i — Vs=o). In order to test in how 
far such a relation can hold we have calculated 
diquark correlation functions in Landau gauge, 

G P , s , c (t) = (D(0)D^t)) . (1) 

Here D(t) denotes one of the states given in the 
second column of Table 1. At large distances the 
correlation functions are expected to decay expo- 
nentially with a mass which is characteristic for 
the given quantum number channel. In Fig. 2 
we show the ratio of correlation functions for 3 C 
diquarks with S = and 1, 

GgosW / G 6i3 0) ~ exp([m 61 3 - 771303] t) . (2) 



Figure 2. Ratio of correlation functions of colour 
anti-triplet states with spin S = and S = 1 
versus Euclidean time t. 

The ratio of correlation function clearly rises ex- 
ponentially at large distances; the slope increases 
with decreasing quark mass. The mass differences 
extracted from fits to this correlation functions 
give, again within simple potential models, the 
difference between the spin-dependent q-q inter- 
actions, 7n 61 3 — 777,303 = (Vs=i — Vs=o). Twice this 
difference is shown in Fig. 1 (open squares). As 
can be seen it is quite compatible, although sys- 
tematically below the results for the N — A mass 
splitting. 

3. OGE VERSUS INSTANTONS 

In order to differentiate between predictions of 
OGE and instanton induced interaction models 
for different quantum number channels one may 
compare the 5 = 1 correlators in 3 C and 6 C chan- 
nels. The ratio G3 16 (t)/G 61 ^(t) is expected to rise 
(fall) with increasing t in the case of OGE (instan- 
ton) induced q-q interactions (Table 1). From the 
results shown in Fig. 3 it is obvious that at least 
for heavy quarks the interaction in the 3 C channel 
leads to lighter excitations, i.e. is more attrac- 
tive. This is in favour of the instanton induced 
interaction models. We note, however, that with 
decreasing quark mass there is a tendency for the 
ratios to flatten. A i-independent ratio in the 
chiral limit thus cannot be ruled out. This could 



3 



G 3l6 (t)/G 6l3 (t) 



i 




0.147 



01 23456789 10 
t 

Figure 3. Ratio of correlation functions of S = 1, 
colour 3 and 6 states versus Euclidean time t. 
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Figure 4. Masses in units of the number of valence 
quarks for different values of k. 



suggest that no bound states exist in these quan- 
tum number channels. 



4. QUARK AND DIQUARK MASSES 

The above discussed analysis suggests that the 
diquark correlator in the 3 C , spin zero channel 
leads to the lightest excitations. The correlation 
function G^(t) does indeed show a rather clean 
exponential decay, which leads to a plateau for 
local masses at distances t ~ 4 for all quark 
masses considered by us. This suggests that 
the diquark is a well localized state. Using the 
string tension to set the scale (y/a = 420MeV) 
we find from an extrapolation to the chiral limit 
m m = 694(22)MeV. 

A similarly stable behaviour is found for local 
masses extracted from the quark propagator. Un- 
like for all diquark masses the plateau is, how- 
ever, reached from below which reflects the non- 
existence of a positive transfer matrix in Landau 
gauge. In the chiral limit we find for the con- 
stituent mass, m q = 342(13)MeV. The mass of 
the lightest diquark obviously is consistent with 
twice the constituent quark mass. There thus 
is no evidence for a deeply bound diquark state 
as it has been found in instanton model calcu- 
lations. A comparison of quark, diquark and 
hadron masses is given in Fig. 4. As can be seen 
1/3 of the nucleon mass lies below the quark mass 



which can be addressed to the contribution of an 
attractive q-q interaction. 

The analysis of diquark correlation functions in 
different quantum number channels shows that 
correlations are stronger in colour anti-triplet 
channels than in colour sextet channels. Nonethe- 
less there is so far no evidence for a deeply bound 
diquark state. One should, however, stress that 
this first exploratory lattice study has been per- 
formed in the quenched approximation on quite 
coarse lattices and still with fairly large quark 
masses. The influence of these approximations 
should be analyzed in more detail in the future. 
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